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Abstract 

In this work, we analyze a Stokes problem arising in the study of the Navier-Stokes flow of a 
liquid jet. The analysis is accomplished by showing that the relevant Stokes operator accounting 
for a free surface gives rise to a sectorial operator which generates an analytic semigroup of 
contractions. Estimates on solutions are established using Fourier methods. The result presented 
is the key ingredient in a local existence and uniqueness proof for solutions of the full nonlinear 
problem. 
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1 Introduction 

In this paper, we are concerned with solutions of the modified Stokes problem 

Av-;uAv + Vg = f on(0,T)xO (1.1) 

V-v = on(0,T)xl7 (1.2) 

v(0,-)=0 onn (1.3) 

S(v,g) = on(0,r)x5ir (1.4) 

D^q\r, = D^q\r„Dlv\r, = Z)tv|r„ for < m < s - 2, 0<k<s-l (1.5) 

for suitable initial data and sufficiently general body forces f . Here Q denotes the set 

n = Dx{0,i), (1.6) 

where £ > and B = {(01,02) £ M? : + 02 < k?} for some radius k > 0. We are primarily 
interested in thin fluid filaments (i.e., where k is small relative to the axial period i) and hence we 
can assume k < 1. Throughout, Cartesian coordinates in M" will be written in the form (ai, . . . , a„). 
Sp denotes the portion of corresponding to the cylinder surface, given by 

Sf = {(01,02,03) G 517 : Oi -I- 02 = K^} , (1.7) 
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To and Ti denote the opposing faces of dQ 

To = D X {0}, = D X {£}, (1.8) 

and s > 2. The quantity v is the (Lagrangian) fluid velocity, q is the (Lagrangian) fluid pressure, 
n is the outward unit normal to O, and > is the (constant) fluid viscosity. The moving 
free-boundary condition is abbreviated by S(v, g) = on Sp, where 

3 3 
S{v,q) = (^qm - n'^{DjVi + DiVj)nj^ _ . (1.9) 
i=i 

Note that the conditions (|1.5p require that solutions be periodic of period i in the as-direction. 
Our objective in this work is to show that this Stokes problem allows unique solutions for given 
initial data and arbitrary T > 0. 

We now briefly motivate how the linear problem (jl.ip - (jl.5p arises in the study of free fluid 
jets. We consider the three-dimensional motion of a jet bounded by an evolving free surface under 
isothermal conditions and without surface tension. The fluid is assumed to be viscous, homoge- 
neous, incompressible, and Newtonian. To model the fluid jet, the three-dimensional incompressible 
Navier-Stokes equations are coupled with periodic boundary conditions in the axial direction as in 
|19j and moving free-surface boundary conditions in the radial direction: 



Dtu + (u • V)u - ^Au + Vp 


= 9^3 


on Q.(t) 


(1.10) 


V-u 


= 


on n{t) 


(1.11) 


{pi - fi{Vu + Vu^)) • N 


= PoN 


on Spit) 


(1.12) 


^3>lr,w=^3>lroW,^3'u|rKt) 


= ^3u|ro{t) 


for0<?n<s-2, 0<k<s-l 


(1.13) 


Dty{t, •) 


= u(t,y(t,-)) 


on Q 


(1.14) 


u(0,-) 


= uo(-) 


on 


(1.15) 


y(o,-) 


= !(•) 


on Q. 


(1.16) 



In this Eulerian description u is the fluid velocity, p is the fluid pressure, y is the fluid parcel 
trajectory map, Pq is the (constant) ambient pressure, and g is the acceleration due to gravity. 
At time t, the fluid domain, free surface, and periodic faces are given by 0,{t) = y(t, 17), SF{t) = 
yitjSp), ro(t) = y(t,ro), and Ti{t) = y(t,T£), respectively. N is the outward unit normal to Q{t) 
and 63 = (0,0,1)^. 

The periodic boundary condition is chosen because it leads to a simpler functional setting and 
avoids all axial boundary layer difficulties while retaining the primary mathematical challenges of 
the problem. In addition, the assumption of periodicity in the axial direction has been successfully 
used to study physical flow phenomena in the numerical simulation of drop dynamics for viscoelastic 
fluid jets in]. 

To obtain a fixed fluid domain, it is useful to shift to a Lagrangian speciflcation of the flow fleld. 
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The problem (|1.10p - (jl.l6p then becomes 



DtVi - fj. Xj,kDk{Xj,mDmVi) + Xi,kDkq = 




for i G {1,2,3} on 


(1 


j,k,m=l k=l 
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on 17 


(1 


j,k=i 








3 

qui- [Xj^kDkVi + Xi^kDkVjjrij = 





tor i G |1, 2, 3| on Di? 


(1 


j,k=i 








DTqW, = D^q\r„Dlv\r, = Dfvlro, I^fxlr, = Dl^\r, 


for < 


m<s-2, 0<A;<s- 


1 (1 




on Vt 




(1 


v(0,-) = uo(-) 


on 




(1 


x(0,-) = 


on Vt. 




(1 



Here v is the Lagrangian fluid velocity, q is the difference between the Lagrangian fluid pressure and 
the ambient pressure, x = y — I is the fluid parcel displacement map, and 5ij denotes the Kronecker 
delta. One consequence of converting the governing equations to the Lagrangian specification is the 
introduction of a priori unknown quantities involving derivatives of the trajectory map y, which 
we denote by Xij{t, a) : — )■ M where 




A = {Xij) = (Vy)-i = D^yi D2y2 1^22/3 ■ (1-24) 



It readily follows from a continuity argument that x ~ for t ^ 1 so that A is approximately equal 
to the 3x3 identity matrix for small times t. Taking Xij = 6ij in (jl.l7p - (jl.23p we obtain, with the 
exception of the initial data, the linearized Stokes problem (jl.ip - (jl.5p . Details of the Lagrangian 
coordinate change for closely related problems are given by Beale [5] and Teramoto |19j . 

To obtain a local- in-time solution for the full nonlinear problem (jl.l7p - (|1.23p . a fixed point 
approach can be used, which requires unique solvability of a slightly more general version of the 
linearized problem discussed here. That general case, however, can be shown to reduce to the 
one treated in this paper. This overarching strategy of studying a modified Stokes problem in 
Lagrangian coordinates and its use in a fixed point argument was championed by Beale in [5] for a 
semi-infinite "ocean" of fluid having a free upper surface and fixed bottom. Teramoto subsequently 
adapted Beale's techniques to gain similar results for a free surface problem involving axisymmetric 
fiow down the exterior of a solid vertical column of sufficiently large radius [19j . 

It is important, however, to note that while the fiuid jet appears similar to the problem con- 
sidered in [19], there are key differences which require that we build upon the work done by Beale 
and Teramoto. For example, unlike the fiuid domains under consideration in [U [U (Sj El [19] , there 
is no stationary surface opposite the free surface to which a Dirichlet boundary condition can be 
assigned. Foremost among the consequences of not having such a condition are the loss of general 
applicability of the Poincare inequality and the loss of invertibility of the modified Stokes operator, 
which is central to the analysis. Moreover, where Teramoto is able to exploit axisymmetry and 
cylindrical coordinates to reduce his problem to two dimensions, the same approach introduces sig- 
nificant challenges in the fluid jet case since the Navier-Stokes equations in cylindrical coordinates 
have singular coefficients when the axis at r = is contained in the fiuid domain. 
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Also in contrast to [31 [H O [6l [19] , we utilize an elegant semigroup approach to the linearized 
problem. This has the benefit of immediately providing a solution to the abstract Cauchy problem 
associated with the problem (jl.l7p - (|1.23p . We improve upon the spectral analysis of the corre- 
sponding operator —A found in [5] and show that it is, in this setting, a sectorial operator which 
generates an analytic semigroup of contractions. Additionally, we are able to establish explicit 
characterizations for both spaces in the modified Helmholtz decomposition of (-L^(r2))^. 

This paper is organized as follows: in Section [2l we introduce the setting of the problem and 
present some preliminary lemmas which adapt and extend standard results from [5] to fit this 
setting; in Section [3l we derive the relevant abstract Cauchy problem and restrict the spectrum 
of the underlying differential linear operator ^4 to a sector in the right half of the plane as well 
as provide estimates on the resolvent operator of A; in Section [31 we establish that —A is the 
infinitesimal generator of an analytic semigroup of contractions and use this solve the linearized 
problem (fLT]) - (fT3|) . 

2 State Spaces and Estimates 

To analyze equations (|l.ip - (|1.5p . we take as our initial fiuid domain (the space occupied by the 
fiuid at t = 0) the infinite cylinder along the aa-axis, 

^oo = {(01,02,03) €M.^ :al+al< k^} , (2.1) 

with free surface d^oo = {(oi, 02, 03) e : of + = k^}. We restrict our attention to flow which 
is periodic in the 03 direction, hence we are interested primarily in functions of the form 

/ = E e hU^^) = <of (^00), (2.2) 

n 

with /„ E H^(J}). In practice however, we will find it more convenient to work with functions over 
a single period. It is natural then to interpret a3-periodic functions on Oqo as being defined on a 
solid torus T C M^. It should be clear that H^{1~) is smoothly isomorphic to the space of functions 
of interest. 

While T is a natural choice for the domain given the periodic setting, we prefer to work in 
the physical space occupied by Ooo- To this end, we notice that there is a C°° diffeomorphism 
from one period of JIqo onto T and consider the bounded set Q = D x (0,£) with boundary 
dil. = Sp U To U r^. While the use of il. in place of iloo does give rise to minor technical issues (as 
opposed to T) concerning the regularity of functions as one approaches the "artificial" corners in 
the boundaries, most of these problems can be dealt with by temporarily exchanging Q for a larger 
subset of f^oo- As such we will occasionally find a use for the set r^i = D x 

Given a spatial domain [/ C and a time interval / C M, the following notational conventions 
are adopted for arbitrary function spaces X{U) and Y{I x U): 



X(C/) = {X{U)f, Y{I xU) = {Y{I X U))\ (2.3) 

X„([/) = {u G X(C/) : V • u = 0} , Y^(/ x [/) = {u G Y(/ x [/) : V • u = 0} , (2.4) 

X = x{n), Y = y((o, T)xn), (2.5) 

°X = {ueX ■.u = on Sf}, °y = {n G y : u = on 5f} . (2.6) 



Here the vector and tensor fields are equipped with the Euclidean and Frobenius norms, respectively. 
To keep the notation simple, we use the following rule: If a function space already has a subscript, its 
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divergence-free subspace will be denoted by appending a a to the existing subscript. Our notation 
is thus closely aligned to the one chosen by Beale in [5]. Spaces not following these conventions 
will be explicitly defined in each instance. We now introduce the spaces fundamental to this text; 
though each assumes as its spatial domain, the extension to is obvious. For the set of functions 
on Q whose as-periodic extensions are continuously differentiable and bounded on we simply 
take = {u\q : u = Er=-oo ^n(ai, a2)e2™'^3/^ G {Tl^}. Note that e C'=(D) necessarily. 
Similarly, we define (or Cp'") to be the set of all such functions which are bounded and smooth 
(or Holder continuous with exponent a) on i^oo- It is clear that the following space is isomorphic 
to H^{T): 

H^ = L= ^in(ai,a2)e2™"^/^ G H^n) : n„ G H''{B) and {u,u)Hk < oo I (2.7) 

I. n=— oo ) 

for A; G Nq, where 

(27rn)^'^ r 

{u,v)Hk= 2^ 2^ £2m-l {Un,Vn)H''~'^iB) and \M = i''^, u) fjk . (2.8) 

n=— oo m=0 

The norms || • Wj^k and || • are equivalent norms on which are actually equal for k G {0, 1}. 
Here || • denotes the standard norm on H'^^Q). It then readily follows that = and 

= {feH^: Dy\r, = Dy\r, for ah < j < - 1} (2.9) 

for A; > 1. We define H^, s G M+, using complex interpolation. Note that, throughout the text, we 
typically use r and s to denote non-integer regularity and k and m when we restrict ourselves to 
integer regularity. 

Instead of the standard Helmholtz decomposition of L^, we take our lead from [5] and pursue 
something slightly different. In particular, to incorporate as-periodicity along with the divergence- 
free condition into the auxiliary space, we choose our decomposition so that can be projected 
onto Hp^. For convenience, we set = Hp^ and introduce the space 

= {v G : Stan(v) = on Sp}, (2.10) 

where Stan = S — (S • n)n is the tangential part of S. Finally, to incorporate regularity with respect 
to time we define the space 

K^{I xn) = H'/'^{I; H^) n H^il; H^). (2.11) 

In contrast to [5], we now provide an explicit characterization of the orthogonal complement 
^pO^± arising in our Helmholtz decomposition of L^. This result will prove important later on. 

Proposition 2.1. The orthogonal complement of P'' in has the characterization (P'^)^ = {Vq : 
q G 'Hi). 

Proof. Let Y = {Vq : q G ^H^}. It is sufficient to show two things: (i) Y is closed in so 
that Y = (Y-*-)"*-, and (ii) P*^ = Y"*-. In order to prove (i), we will first need to show that the 

orthogonal complement of X = '^C^ in has the characterization X"*- = {Vq : q G H^}- Let 
q G Hp,u G X. There exist G '^C^ such that Ufc — )• u in L^. Integration by parts yields 

(Vg,u)L2= lim (Vg,Ufc)L2 = lim / gu^ • es + / qu^ ■ {-e-i) = 0. (2.12) 
fc-5.oo fc^oo Jto 
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Thus Vg G X-*-. Conversely, let w G X-*-. Then, in particular, (w, u)l2 = for all u G C^. Thus 
there exists p G such that w = Vp, see |18l pp. 10-11]. Now consider 

u=( )g°C-, (2.13) 

\u{ai,a2)/ 

where u G C^(D) is arbitrary. Then, applying integration by parts, we obtain 



= (Vp, u)l2 


(2 


= pu e3+ pu- (-ea) 


(2 






= pu — pu 


(2 






= ip\r, -p|ro,'«)L2m) • 


(2 



Since u is an arbitrary element of a dense subset of L^(D) (see \12\ p. 13]), this implies that 
pWi = pIto on L^(D). Hence p G 



With this characterization in hand, we can now prove (i). Let Qk G such that Vg^ — t- f G L^. 
Since Vq^ G X^, we have 

(f,u)L2 = lim (V(7fc,u)L2 =0 (2.18) 
for all u G X. Hence f G X"*- and so there exists p G Hp such that f = Vp. Notice that for n 7^ 

||(^fc)n -Pn 1 1^1(15) ^^^X] ^("7") ll(^fc)n-Pn|lHl(IO) +5^ll^i((^fe)n-Pn)|li2(e)j (2.19) 

<i^\\Viqk-p)\\l2. (2.20) 

Thus (gfc)„, — )• in -ff"'^(ID'). Since {qk)n G -^^0(1^)5 a closed subspace of H^{D), we obtain p„, G -ffp (B) 
for n 7^ 0. For n = 0, applying the standard Poincare inequality yields a constant C > such that 

||(4)0 - (grn,)o||^i(o) < C\\V{{qk)o " (gm)o) ||?i2(iD,))2 < Ce^\\V{qk " 9^)11^ (2-21) 

which implies that {qk)o converges in Hq{0). Moreover, the limit is necessarily po + A, for some 
A G M, since it is readily seen that {qk)o converges to this limit in the weaker L^-norm. Thus f = Vq 
where g = p + A G ^H^. Hence Y is closed in L^. 

Finally, we show (ii). Let u G and 93 G C^. Then 

= (Vip, u)l2 = - I ifiV^). (2.22) 
Jn 

Hence V • u acts as a bounded linear functional on and can be extended to all of by density. 
This unique operator must be the zero functional and thus u G P'^. Conversely, let v G P*^. Since 
= Y © Y^, there are q G ^H^ and v G Y^ such that v = v + Vg. Taking the divergence of 
both sides of this equation yields Aq = and, by Lax-Milgram, q must be the unique solution of 
this equation in ^Hp. Thus q = and v = v G Y-*-. Thus P*^ = Y^ and the claim follows. 

□ 
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The following are analogous to results in [5] and are readily shown using similar techniques. 
They are given here explicitly for the reader's covenience. 

Proposition 2.2. Let P be the orthogonal projection of onto P''. 

(1) For s > 0, we have -PHp = P^ and -PIh^ : Hp — P^ is bounded. 

(2) -P|kj : Kp — )■ Kp is bounded with norm bounded independent of T. 

(3) Suppose s > 1. If / G H^, then there is a unique / G such that 

P{Vf) = Vf, f\s, = f\s,, and A/ = 0. (2.23) 

We will see that, just as in [5], many crucial quantities can be cast as solutions of a particular 
problem involving Laplace's equation. Adapting the boundary conditions to reflect periodicity and 
the absence of a fixed bottom surface, the relevant problem in our setting takes the form 

An = / in 0, n = on Sp, D^u\rt = -Cfwlro for k G {0, 1}, (2.24) 

where / G H^~^ is given. The following result demonstrates that (12.24p has a unique solution and 
provides an estimate for it in terms of the inhomogeneity /. We note that the provided proof does 
not draw from the corresponding proof in [5]. 

Proposition 2.3. For / G H^~^,s > 2, there is a unique solution u G of Au = f on Q. 
Additionally, there exists C > 0, independent of /, such that ||ti||_f/s < C'||/||j^s-2. 

Proof. Let / = /nC^'^*"'^^/^. We first consider the boundary-value problem, L^n = —fn on B 
with n = on 9D, where L„ and its associated sesquilinear form (S„ : Hq(I})) x Hq(B)) — )■ C) are 
given by 

LnU = -An+ ( — J u (2.25) 

/27rn\ ^ 

Bn[u,v\ = (Vn, Vu)l2(io,) + —— (u,u)i2m). (2.26) 



Clearly, Bn is continuous and coercive on Hl{B), thus we can apply Lax-Milgram to obtain a 
unique weak solution, Un G Hq{J}). The construction u = then our candidate for 

the solution of the boundary- value problem in O. We now restrict our discussion to the case when 
s = A; G Z. Given the regularity of we can immediately conclude that each Un G H^{B)) is a 
strong solution. Our goal is to show that u G H^. First we obtain some preliminary estimates for 
Un where n 7^ 0: 



Bn[Un,Un] = {-fn,Un)L^m (2.27) 

/ 27rn\ ^ 

L2(D) 



|V''^n|lL2(D) + ~~r- ll^nllL2(D) - ll/nllL2(ID))ll^^n||L2(p) (2.28) 



2 



||'"n||L2(p) < WfnWL^iO) (2.29) 

Notice that from this estimate we can conclude 
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This gives us an estimate on the lowest order terms in the Hp-norm. For the highest order terms, 
standard elhptic regularity theory (e.g., see [9], p. 323) provides an estimate of the form 

ll''in||j?''-(D) < C'l||/n||_H'*-2(I[I))) (2.31) 

though the constant Ci here generally depends on the coefficients (and hence n) of -L„. However, 
upon closer inspection of the proof of this result (e.g., in [9]), we observe that our above estimates 
on ||iin||L2(iD)) can be used in place of the usual L°°-estimate on the coefficient (27rn/£)^ of L„. In 
consequence this allows Ci to be chosen independently of n. Therefore 

E < C! < < oo. (2.32) 

n n 

Finally, we must show that the intermediate order terms in the i/p-norm are also summable. 
Exploiting complex interpolation between H^(J])) and H'^iJ}) and Young's inequality, we obtain for 
each < m < /c 

[['""llgfc-'^m) - / . \~f j (^ll^"llL2{D)ll'"'"ll/ffc{D) j (2.33) 

/r, \ 2m(fc-2)/fc 



k/m 



+ - 



k — m 



k 



II ? ||2(fc— m)/fc 
\\jn\\jjk-2(o) 



m ( 27rn 



2(fe"2) 



ll/nl 



L2f 



+ 



k — m , 



k 



k/{k—m) 



<Cs\\f\\l, 

Thus u G with ||n||2 < Ci{k + 1 
Interpolation then provides the remaining cases. 



(2.35) 

(2.36) 

(2.37) 
(2.38) 



, , which completes the proof for integer values of s. 

□ 



3 A Spectral Result 

Our first goal is to use the modified Helmholtz projection P to rewrite the problem (jl.ip - (jl.5p in a 
variational form which has the velocity as its only unknown. First we notice that for any solution 
(v, q) of the problem, ()1.2p implies v(t) G for each t. Since it is readily seen that P commutes 
with Dt, applying P to (jl.ip yields 

Av - AiPAv + Vgi = Pf , (3.1) 

where Vqi = PVq (with Ag'i = on ^2 and qi = q on Sp) by Proposition 12.2( 3) . This application 
of P removes the indeterminacy of the pressure term in the sense that the value of qi is uniquely 
determined by v. 
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Lemma 3.1. Suppose s > 2 and (v, g) G Hp x ^ satisfies ()1.4p . Then there exists a bounded 
Unear operator Q : Hp — t- -ffp"^ mapping v i— t- gi where gi is the function provided by Proposition 
with Vqi = PVq. 



Proof. To see this, we use the fact that q and qi agree on the free surface and observe that (II. 4p 
imphes that the normal component of S(v, q) must vanish on Sp. Hence S(v, gi) • n = implies 
qi = 2/iK"^ Si j=i CLiajDjVi on Sp- Given v G Hp, we note that / = 2/xk~^ Yld j=i o-iO-jDjVi € 
ifp~^. For s = 2, we can apply Lax-Milgram to obtain the existence of a unique weak solution 
qi G Hp of the problem 

Aqi =0 on O (3.2) 

qi = f on S'f (3.3) 

Dlq,\r, = Dfgilro for fc E {0, 1} (3.4) 

with ll^ill//! < Ci||v||h2 where Ci > is independent of v. For s > 3, we consider the problem 

A(p = -A/ on n (3.5) 

(/> = on Sf (3.6) 

DgVir, = I^sVlro for k e {0, 1} (3.7) 



3 
P 



which has a unique solution cp G ^H^ ^, by Proposition 12.31 satisfying < C2||A/||^, 

for some C2 > which is independent of v. Finally, we set qi = (p + f ^ H^^^ and observe 
that < C2||A/||^s-3 + ll/llji^s-i < CsIIvIIh^- Interpolation now yields the claim for the 

remaining values of s. It readily follows that the constructed operator is linear in v. □ 

We now take the general approach used in semigroup theory by treating (jl.ip as an abstract 
ordinary differential equation with respect to time whose solution is, for each value of i, an element 
of the appropriate function space (V'^) on If we define an operator A : V'^ — >■ p^"^ by 

= -fiPAv + VQv, (3.8) 

the problem (jl.ip - (jl.5p takes on the form 

v + Av = Pf on(0,r)xfi (3.9) 

v(0,-)=0 on (3.10) 

Notice that (|1.4p is satisfied since our construction of Q ensures that the normal component of 
S(v, g) will vanish on Sp. The operator A is a modification of the standard Stokes operator. 

We tackle the matter of determining the spectrum of —A first. Unfortunately, in contrast to 
the problems treated in O HI [5l [H [19], A is not injective with our boundary conditions (implying 
that lies in the spectrum of A) since ^(v + c) = Av for any constant vector c. This, combined 
with the inability to apply the Poincare inequality in general, makes the problem of determining 
the spectrum more challenging here than in the aforementioned cases. Restricting the spectrum of 
A to a sector in the right half of the plane and providing estimates on the resolvent operator, the 
following theorem is a key result of this work. 

We will show that p{A) contains all A € C such that |Im(A)| > Re(A). Given g G ps~2 ^^^^ ^ 
with |Im(A)| > Re(A), we find a unique solution {\r,q) G V* x H^"^ of the problem 

-//Av-Av + Vg = g (3.11) 
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along with (jl.2p . (jl.4p . and (jl.5p . To see that this is equivalent to the statement about p{A), 
suppose that there exists v € such that [A — A/)v = g. Using our decomposition of L^, 
we find qo G such that Vqq = — P)Av. Setting q = Qv + goj we obtain (j3.1ip . It is 
now straightforward to verify that (v, q^) also satisfies (jl.2p . (jl.4p . and (II. 5p . Conversely, given a 
solution (v, q) of the stationary problem we can apply P to (j3.1ip to obtain {A — AI)v = g. Hence 
{A — A/)v = g has a unique solution v if and only if the problem p.2p . (|1.4p . (II. 5|) . (I3.11|) has a 
unique solution (v, q). 

Theorem 3.2. For g G P''"^, the problem (fOj) . (fL4l) . (fT3]) . (fXTTI) has a unique weak solution 
(v,g) G X for each A G C such that |Im(A)| > Re(A). 

Proof. Notice that for any v G P^, v satisfies (jl.2p and the portion of (jl.Sp referring to the 
velocity. The free surface condition (11. 4p is not necessarily satisfied though and will need to be 
incorporated into a variational formulation directly. To that end, let us consider the sesquilinear 
form (•, •) : P^ X P^ — ^ C defined by 

3 „ 

(v,u) = -A(v,u)l2 + ^^2 (^i^^ + DiVj){DjUi + DiUj). (3.12) 

Now suppose u G Hp, v G P^, (7 G and observe that 

/ (-/iAv - Av + Vg) • u = f / V Ai;iIIi ) - A(v, u)l2 + / g(u-n)- /gV-u (3.13) 
Jn \Jn . J Jdn Jn 

= 11^2 DjUiDjVi + / Uiiqui - /i DjViUj) 
^Jn ^Jan ^ 

-A(v,u)l2 - [ qV-u (3.14) 
Jn 

= (v,u)+ / S(v,g) • U + /i / UiDiVjHj 

Jdn ~i Jan 

- / DjUiDiVj- qV-u (3.15) 
Jn Jn 

= (v,u)+/" S(v,g) -u- /" gV-u + ^V /'lI,A(V-v) (3.16) 
JOf^ Jn Jt7 



(v,u)+ / S(v,(?)-u- / gV-u. 



(3.17) 



Notice that the pair (v, q) currently satisfies ()1.2p and (jl.Sp . If we suppose that (v, q) additionally 
satisfies ()3.1ip and (jl.4p . then we obtain 

(g,u)L2 = (v,u)+ / S(v,g)-u+/ S(v,g)-u=(v,u) (3.18) 

for all u G P^. Thus (v, u) = (g, u)l2 can be seen as a weak formulation of the full problem which 
does not involve q. In an effort to apply Lax-Milgram, we verify that the sesquilinear form is both 
continuous and coercive. Applying Holder, |(v, u)| < C||v||hi ||u||hi , where C > depends on 
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/i and A. Hence the sesquilinear form is continuous. That it is also coercive follows from Korn's 
inequality: 

Kv,v)|2= JjD,v, + D,v,\'-Re{X)\\v\\l)j + {lm{X)Ml)' (3-19) 

DjVi + D,vj\^ - Re(A)||v||22 + |Im(A)| • \\v\\l}j (3.20) 

(3.21) 

For Re(A) > 0, line (ICTD imphes 

Kv,v)|2>Im(A)2||v||^. >i|An|v||^.. (3.22) 

Moreover, the same estimate can be obtained for Re(A) < since the right-hand side of line (j3.19p 
then expands to something of the form (j) + |Ap||v||^2 where (p > 0. Since the sesquilinear form 
satisfies the conditions of Lax-Milgram, we obtain a unique weak solution v G of (USD, P^ . 
(jl.Sp . (jS.lip along with the estimate 

I|v||hi < C\\gU2. (3.23) 

We now seek an associated pressure, q, of v. Recall that an associated pressure need only satisfy 
()3.1ip in the sense of distributions (i.e., when tested against arbitrary u G CJf ). As with the 
velocity, we begin by finding a weak formulation for the pressure. Notice that, for q £ with 
(v, g) satisfying (jl.4p . we obtain from (j3.17p that 

/ (7V • u = (v, u) - (g, u)l2 (3.24) 

for all u G Hp. Using continuity of the sesquilinear form we obtain immediately that the right-hand 
side is a bounded linear functional in u, F : — t- C, which vanishes when V • u = 0. From [15] . 
we know there is a unique q £ L'^ such that 

F = Vq and q = 0. (3.25) 



It is now straightforward to verify that q = —q satisfies (j3.1ip in the distributional sense and hence 
is an associated pressure of v. It is uniquely determined under the additional condition j^q = 0, 
but otherwise is unique only up to a constant. □ 

The proof of the preceding theorem allows us to draw the following conclusion about the case 
A = 0. 

Corollary 3.3. For g = 0, {v,q) G x is a weak solution of the problem (fL2]l . (fLi|l . (fT3]l . 
(j3.1ip with A = if and only if v and q are constant. 

Proof. Reasoning as in ()3.19p . we have 

= Kv,v)| = ^ V / \DjV^ + DiVj\^ > 2/i V / \D,Vi\\ (3.26) 
Hence v and a fortiori q are constant. □ 
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Theorem 3.4. Let s > 2. Then a{A) C {A G C : |Im(A)| < Re(A)}. Moreover, for A with 
|h7i(A)| > Re(A) and |A| > e > the resolvent operator R{X; A) = {A — XI)^^ satisfies 

\\R{\;AMns < C{\\g\\^s-. + (1 + .-^)(|A| + (3.27) 

for all g G P^~^. Here C > is a constant which is independent of A, e, and g. 

Proof. We will now demonstrate that the weak solution provided by Theorem 13.21 can, in fact, be 
made into a strong solution of the problem (jl.2p . (jl.4p . (jl.Sp . (j3.1ip . In order to avoid the lack 
of regularity due to the "artificial" corners in our domain, we turn to the equivalent problem of 
finding a weak solution (vi,qi) of the problem (jl.2p . (|1.4p . (jl.Sp . (j3.1ip on the larger domain Jli. 
By choosing qi such that J^^ gi = we can ensure that (vi, qi) is simply the periodic extension of 
(v, q) to Oi. Then it follows from standard results [ini that (vi, gi) has the additional regularity we 
seek on compactly contained subsets of Oi. To obtain regularity all the way up to the boundary, we 
follow the approach in [16] which is applicable to the boundary provided that it is smooth locally. 
Thus the pair (vi, gi) has the desired regularity near 5ir up to and including the intersections with 
Tq, Ti since these regions occur on a smooth portion of the free surface on Oi. It follows that 
G Hf^c(^oo) and qp E iJ/^^(Ooo), hence v G and g E iJj^. 

To see that (v, q) provides us with a strong solution of our problem, we only need to verify that 
()1.4p and (j3.1ip are satisfied. Using (j3.17p . for all u G we have 

(-^fAv- Av + V<?-g,u)L2 = / S(v,(7)-u= / S(v,,?)-u. (3.28) 

Jan Jsf 



Taking u G "^C^ implies that — ^uAv — Av + Vg — g lies in the orthogonal complement of *^C^ ^ . 
It was demonstrated in the proof of Proposition 12.11 that this orthogonal complement consists of 
the gradients of functions in H^, so that —fiAv — Av + Vg — g = Vp for some p G H^. However, 
(j3.17p now yields 

(g, u)l2 = (-/iAv - Av + V{q - p), u)l2 = (v, u) + / S(v, q-p).u- [ {q-p)V-u (3.29) 

Jan Jn 

for all u G Hp. Restricting u to and exploiting (13.25P reduces this to f^pV -u = 0. Integrating 
by parts, we see that Vp ■ u vanishes for arbitrary u G CJf . Since this is a dense subset of L^, 
Vp = and q satisfies (|3.1ip . All that remains is to show that (jl.4p . the free surface boundary 
condition, is also satisfied. From (j3.28p we now immediately obtain 



/ S(v,g)-u = (3.30) 

Jsf 



for all u G P^. Following the lead of [TTj, we localize to a neighborhood S C Sp and construct 
u G P^ such that u\sp = (S(v, g) — (S(v, g) • n)n)(/) where (p is a smooth nonnegative function 
vanishing outside S. Then 

/ S(v, <7) • u = / |S(v, q) - (S(v, q) ■ n)n|2^ + (S(v, q) ■ n)n • (S(v, g) - (S(v, g) • n)n)<A (3.31) 
JSp Jt, 

= [ |Stan(v)P^ (3.32) 

= (3.33) 
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implies that Stan(v) = on S. Since S was chosen arbitrarily, we obtain S{v,q) = {S{v,q) ■ n)n 
on Sp- Let 6{v, q) = q — 2/iK~^ ^ aiajDjVi G H^. Since 0{v, q)\sF = S(v, q) ■ n, (I3.28|) yields 

/ 9{v,q)n-u = [ V9{v,q)-u = (3.34) 

Jan Jn 

for all u G P^. By density, V0(v, g) G (P*^)"*" and 0{-v,q) = p + ui for some p G "^-ffp and w G M. 
Since this implies S(v, q) ■ n = q — 2/ik~^ Y1 cuajDjVi = a; on Sp, we take q* = q — co and obtain a 
unique strong solution (v,g*) G x i^^^ of the problem (fO|) . pTi]) . (fT3|) . ([XTT]) . 

To further increase regularity, we turn to the standard a priori estimates of Agmon, Doughs, 
and Nirenberg (ADN) [2]. Here it is useful to work on T rather than $7. The associated problem is 
readily seen to be uniformly elliptic in the sense of ADN with boundary conditions satisfying the 
complementing condition. The a priori estimates in T then lead to the following estimate in 0: 

3 3 

ll^*lli^r' +^"''^'"-f^p ^'^^J2Mh;-^ (3-35) 

for a positive constant Cx which depends on A. Thus (v, g*) G V x H^~^ is the unique solution 
of (^2^, (mi), ([LSD, <|3TT]) and a{A) C {A G C : Re(A) > 0}. Now all that remains is to show that 
the resolvent estimate (j3.27p is satisfied. From (j3.35p we obtain the estimate 

Mi,s< (Ell^'^^-ll^pj <3C|||g||^._. =3C7|||(A-A/)v||^._.. (3.36) 



Thus we have 



|v||h5 <ci||(^ + I)v||h.-2 (3.37) 

< ci (||(^ - A/)v||h.~2 + (|A| + 1)||v||h.~2) (3.38) 

< C2 (llgllHr^ + (|A| + l)||v||{^7'^/^||v||J^//) , (3.39) 



where ci and C2 are positive constants which do not depend on A. Here we have used complex 
interpolation between and Hp. Finally, we apply Holder to (I3.22P which yields 

IIsIIl^ > ^I|v||l2. (3.40) 

Now let us restrict ourselves to |A| > e for arbitrary e > 0. If s = 2, then (j3.39p and (13.40p yield 
(j3.27p directly. Otherwise, we can apply Young's inequality to (j3.39p obtain 

||v||h5 < C2 f||g||H-2 + C3(|A| + 1)*/2||v||l2 + ^I|v||h^) (3.41) 



<C4(||g||H-2+(|A| + l)^/2||v||L2j (3.42) 

< C5 (hWui-^ + (1 + e-'m + l)(^-2)/'||g||L0 , (3.43) 



where C3, C4, and C5 are positive constants which do not depend on A and e. Since v = R{\;A)g, 
this completes the proof. □ 

Since Hp is compactly embedded in Hp~^ for s > 2, Riesz-Schauder theory and Corollary 13.31 
imply the following result. It follows that the kernel of A contains constants only. 

Corollary 3.5. cr{A) consists of isolated eigenvalues of finite multiplicity. Moreover, the eigenvalue 
of ^ has multiplicity 1. 
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4 The Inhomogeneous Cauchy Problem 



We can now show that —A is the infinitesimal generator of an analytic semigroup of contractions. 
This is the main step involved in constructing solutions to the linear problem (jl.ip - (jl.5p . We refer 
the reader to [8] for standard results in semigroup theory. 

Theorem 4.1. The operator — A, with domain V^, generates an analytic semigroup of contractions, 
J{t), on po with ||J(t)|| = 1. 

Proof. As we seek to apply Lumer-Phillips, we begin by showing that —A is dissipative. To do 
this, we must improve (slightly) upon the estimate provided by (j3.40p . For A < 0, we obtain 

3 „ 

l(g,v)L2| = Kv,v)|=-A||v||2, + ^ ^^ J \D,v, + D,v,\'>-XMl,. (4.1) 

Dissipativity now follows using the Holder inequality. Since ^ + / is surjective by Theorem 13.41 and 
P'^ is reflexive (as a Hilbert space), we can apply Lumer-Phillips to obtain that is dense in P'' 
and — A generates a Co semigroup of contractions, J{t), on P''. As the generator of a Cq semigroup 
of contractions, —A is closed (see Theorem II. 1.4 in [8], for example) and using (j3.40p together with 
Theorem 12.31 from [13] we see that J{t) is actually an analytic semigroup on P". Now, since J{t) 
is a semigroup of contractions, we have || J(i)|| < 1. However, is contained in the point spectrum 
of —A (see the discussion preceding Theorem 13. 4p . Hence || J(t)|| = 1 as required. □ 

With this semigroup result in hand, we are finally ready to solve the inhomogeneous linear 
problem (jl.ip - ljl.Sp . Theorem 14.11 immediately provides a solution to the Cauchy problem ()3.9p - 
(j3.10p and makes the Paley- Wiener theory utilized in [5] unnecessary. Here we abbreviate the sets 
G = {0,T)xn and dGp = (0, T) x Sp. 

Theorem 4.2. Let 3 < s < 4, T > 0, and f G K^-^ such that Pf (0, •) = 0. Then the problem 

(|l.ip - (jl.5p has a unique solution {^r,q) such that v E Kp, Vq G Kp~^, and q\sp G Kp ^^"^{dGp)- 
Moreover, this solution satisfies 

I|v||k^ + ||Vg||K-2 + Msf\\k;-^/\qgp) ^ ^llfllKr^' (4-2) 
where C is a positive constant which is independent of T and f . 

Proof. First we notice that Pf € C'^'^''~'^^/^([0, T]; P^) by the Sobolev Embedding Theorem. Com- 
bining Corollary 4.3.3 and Theorem 4.3.5(iii) from |13j . the abstract Cauchy problem 

V + Av = Pf (4.3) 
v(0, •) = (4.4) 

has a unique strong solution v G C^'^*~'^^/^([0, T]; P*^), with v(t) G for each t G [0, T]. Here we 
are exploiting the fact that —A is the generator of an analytic semigroup on P''. Note that v is a 
strong solution in the sense of semigroups, i.e. , v is differentiable almost everywhere on [0, T], with 
V G L^((0,r);P°), such that v(0, •) = and v(t) = -Av(t) + Pf{t) almost everywhere on [0,r]. 
In fact, V is a classical solution in the semigroup sense since it is continuously differentiable with 
respect to time. 

To show that v G Kp, we reconsider the abstract Cauchy problem (now with a new unknown 
variable v) from another perspective. We begin by applying the periodic analog of Lemma 2.2 from 
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[5] in order to extend Pf to Kp~^(R x Q) in such a way that the extension is bounded independent 
of T and vanishes for t < 0. Multiplying through the abstract Cauchy problem by the weight 
w{t) = e~* and taking Fourier transforms in t, we obtain 

Tum) = + (1 + im^'^upm)- (4.5) 

Since it is clear that J^^(Pf)(^) G P'^~^, this uniquely defines -F^(v)(^) e V* by Theorem 13.41 
Making use of the Fourier transform characterization of '^-spaces for s G M"*^ (e-g-, see [1]) and 
the fact that Fourier transforms are unitary transformations, we have 

l|v|lK|^(Exn) ^ 2 (^||v||l2(]r.^s) + I|v||hV2(R;L2)) (4.6) 

= 2 + mhiu;Hs) + 11(1 + er^'^um + mi^^;L^)) (4.7) 
= 2 [ (\\Tum + mh^ + {i+er^'\\^um + mh) (4.8) 

Applying the resolvent estimate (|3.27|) to the first term of the integral, we obtain 

w^^um + i)iIh^ < ci [w^^ipm + i)iiHr + 2(11 + ^(e + i)i + i)^^~'^/'\\T^{pm + i)iil2 

(4.9) 



< C2 [WT^pm + i)iiHr^ + (v^rTTe + iF+ i)'-'ii-^«,(^f)(e + i)iiL^ 

^(4.10) 

< C2 (ii-F^(pf)(e + 1)11^.-2 + (371+^2)"' WTUPm + mh'^ (4.11) 

< C3 (||-F^(Pf)(e + 1)11^.-2 + (1 + ^2){-2)/2 ||^^(pf)(^ + i)||2^) (4,12) 

where ci, C2, and C3 are positive constants which are independent of ^ and f. Similarly, we can 
apply estimate (|3.40p to the second term of the integral to get 

(1 + er^^^um + mh < 2(1 + er^'^ii + + ir'n^^ipm + mh (4.13) 

Combining these estimates yields 

iivii^.(Mxf^) < C4 / \\TUpm+mns-. + {i+ef~'^^'\\:Pn.ipm+mhdc (4.16) 

= C4 (||-F.(Pf)(^ + l)ll^2(j,^^.-.) + 11(1 + e)^'-'^^'j'n,{Pm + l)ll^.(M;L2)) (4.17) 
= C4 (ll-ff|lL2(R;H»-2) + ll-f'f|lH(-2)/2(R;L2)) (4-18) 

^ ^4||Pf||^.-2(^,^), (4.19) 

where C4 > is a constant which is independent of ^ and f . By uniqueness, we must have v = 
v\g G Kp. We now seek a suitable q so that (v, g) is the unique solution of (jl.ip ^ ijl.Sp . For fixed 
t, this amounts to finding a unique q G H^~^ such that 

Vq = ^iAv + + f - Pf on (4.20) 

q = Qv on Sf- (4.21) 
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Since s > 3, this is easily accomplished by taking the divergence of the first equation and applying 
Proposition 12.31 All that remains is to show that (v, q) satisfies (I4.2p . To estimate q we first notice 
that 

Vq = - P)Av + VQv + (/ - P)f . (4.22) 

The only term which we do not yet know how to estimate is VQv. However, since AQv = on 
and Qv = (p on Sp where cp = 2^/^^^ Xli j=i ^i^j-^i^i ^ Hp~^^ if follows from Proposition 12.2( 3) 
that VQv = P{V(t)). Then by Proposition [22];2), 

l|vgv|iK.-2 = ||p(v,^)||k-2 < c5||v,^||k.-2 < csHvUkj, (4.23) 

where C5 and cg are positive constants. Similarly, since Q was constructed so that q = Qv on Sp, 
11^-3/2(3^^) = IIQv|s^||^.-3/2(g^^^ < C7||Qv||^.-i < csIIvIIk^, (4.24) 
where C7 and cs are positive constants. Thus, combining estimates, we obtain 

I|v||kj + \\Vq\\j^s^-2 + \\q\sAK;-y\dGp) - (h^Hk^ + W^Wk^-'') (^-2^) 

< eg (||v||K-(]Rxn) + l|f IIk;-^) (4-26) 

< cio (l|h||K-2(Rxn) + l|f|lKrO ^^•^'^^ 

< cii||f||j^.-2 (4.28) 

where cg, cio, and cn are positive constants which do not depend on f (or T). □ 

5 Concluding Remarks 

With Theorem 14.21 in hand, it is now straightforward to follow through the fixed point approach 
outlined in [Sj, with minor revisions, to obtain the following local existence result. For details about 
these modifications we refer the reader to [7j. 

Theorem 5.1. Suppose 3 < s < |. For any uq G V*~^ there exists T > 0, depending on 
||uo||jjs-i, so that the problem (|1.17|) - (|1.23|) has a solution (v, g) with v G Kp,g G Kp '^^'^{dGp), 
and Vq G K^-^. 

Since the same arguments can be successfully applied to yield a similar local existence result 
when the initial displacement is taken to be nonzero in (|1.23p . uniqueness of solutions can then be 
shown to follow in the standard way. It is also a simple matter to prove that the unique solution 
given by Theorem 15.11 is axisymmetric provided that uq is taken to be axisymmetric. [7] contains 
additional details. 
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